Crossover from BCS superconductivity to BEC of pairs: The role of the lifetime of 

the pairs. 
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The understanding of an electron gas with short coherence length pairs formed by an attractive 
interaction is believed to be one of the major keys to our theoretical knowledge of the high-T c - 
superconductors. Mainly the deviations of the cuprates from usual metallic Fermi liquid behaviour 
already in the normal state like e.g. a linear resistivity or the observation of a pseudo gap can result 
from electron-electron correlations. 

We therefore investigate the negative U Hubbard model in two dimensions at low densities using 
the T-matrix approximation. In the non selfconsistent formulation of the theory the system always 
shows an instability towards Bose condensation of pairs into an infinite lifetime two-particle bound 
state. If the calculations are performed selfconsistently pair-pair scattering is included which causes 
the pairs to have finite lifetime. The physics of these finite lifetime pairs is discussed. 
Keywords: negative-U Hubbard model, two-particle bound states, pseudogap, non Fermi-liquid 
properties 
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I. INTRODUCTION: 

More than a decade ago the high-T c superconductors 
have been discovered. Even if many features of these 
materials like e.g. the mainly two dimensional transport 
in the hole doped CuC-2 planes are understood the mi- 
croscopic mechanisms and models which lead to a wider 
understanding still need to be investigated. 

Already in the normal state above T c there are two im- 
portant issues which need to be explained. The first one 
(I) is the linear resistivity. In a normal metal the low- 
temperature dependence of the resistivity is governed by 
either impurity scattering which gives a constant term or 
by a Bloch-phonon contribution (~ T 5 ) or in very clean 
Fermi liquids by the Landau electron-electron scattering 
of electrons at the Fermi surface (~ T 2 ). Only at high 
temperatures above the Debye temperature which is of 
the order of 300-400 K normal metals show a linear re- 
sistivity due to phonon contributions. 
In the high-T c 's however the optimal doped samples show 
a linear resistivity over the full temperature range, al- 
ready well below the Debye temperature 0. 
The second feature (II) in the normal state which needs 
to be explained is the occurrence of a pseudo gap. It has 
been found in many experiments (first in NMR 0) and 
also with other methods e.g. tunneling experiments 
This pseudo gap is found in the underdoped materials 
were already at temperatures above the superconduct- 
ing transition temperature T c a gap opens at the Fermi- 
surface. 

Also in the superconducting state there is besides the 
high transition temperature another special feature of 
the cuprates. This is the short coherence length of the 



pairs. Usual weak coupling BCS superconductors have 
coherence lengths £ of several hundred lattice constants 
(~ 1000 A) which go together with high quasiparticle 
densities. Therefore many Cooper pairs (~ 10 6 ) over- 
lap and the two essential conditions for superconductiv- 
ity, pairing and phase-coherence, occur simultaneously 
at the same temperature, the mean-field (BCS) T c . 
In the high-T c superconductors however the coherence 
length is very small (£ w 10 A) of the order of only 3-4 
lattice constants. Further the quasiparticle densities are 
very low since the materials are close to an antiferromag- 
netic insulator to metal transition. Therefor the pairs 
barely overlap and might even be treated as spatially 
well separated pairs. This is believed to be the reason 
H why pairing (opening of a pseudo gap?) and phase 
coherence (T c ) might occur at different temperatures. 
All these features above motivate the study of a model 
system which allows to investigate the transition between 
usual BCS superconductivity to a situation were tightly 
bound pairs occure. The simplest model system which 
can describe this is the Hubbard model with an attrac- 
tive interaction between the electrons. 



II. MODEL: 



The negative U Hubbard model is 



(1) 



<hJ> 



were the kinetic energy is given by the transfer term t 
and the Hubbard attraction U is chosen to be negative. 
In the weak coupling case in three dimensions the Fouric 
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transform of eq. ([j]) describes just the starting Hamilto- 
nian for the BCS theory. 

In the zero density and strong coupling (atomic) limit 
(n/U = 0) the Hamiltonian is solved by local pairs of elec- 
trons whose energy is lowered by the attractive energy U. 
Such pairs of electrons can be considered as Bosons. 
In the intermediate coupling regime (|C/| ~ bandwidth 
W) and for low densities n short coherence length pairs 
which barely overlap can be described. This is the pa- 
rameter regime we are interested in in the current work. 



which determines T c in the BCS theory. 
In the second case (II) the chemical potential is below 
the one-particle continuum. When considering \° this 
can always lead to a bound state at K = in ID and 
2D. In 3D this condition only leads to a bound state 
below the one-particle continuum if \U\ > U cr u. Such a 
bound state can be populated by pairs of electrons which 
obey Bose statistics and therefore Bose condensation can 
occur. A Fermi surface will be lost at low temperatures. 
The K-dispcrsion of such a bound state is shown in fig. 
g (a). 



III. T-MATRIX: 

An approximation which allows to access the low den- 
sity regime and which also describes pair formation is the 
T-matrix approximation. It is often also called ladder- 
approximation or Bruckner-Hartree-Fock theory. In this 
theory the infinite sum over all scattering events between 
two electrons enters in the vertex function V . 



f(K,ifi n ) 



U 



l-f/x(K,ifi r 



(2) 



With K being the total momentum of a pair and 
VL n Bosonic Matsubara frequencies. The susceptibility 
X(K, iO n ) is given by the product of two one-particle 
Green's functions G(k, iu) m ). 

x(K, iO n ) = — ^ G ( K _ k ' ~ iw m )G(k, iuj m ) 



m,k 



(3) 



Spin indices are left out since only singlet pairs can be 
formed with an on-site attraction. 



A. Thouless-criterion: 

The physics which is described with the T-matrix ap- 
proximation depends crucially on the appearance of a two 
particle bound state which is given by a zero of the de- 
nominator of eq. (^|) reaches the chemical potential first 
at K = and is called the Thouless instability M. 
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(4) 



Here we have to distinguish between two different cases: 
In the first one (I) the chemical potential is in the one 
particle continuum. One has a Fermi surface. This hap- 
pens in the weak coupling 3D case when a BCS instability 
is formed at the Fermi surface. The Thouless instability 
occurs when the system is cooled down towards T c . Eq. 
(||) is (if x — X° is build up from noninteracting one- 
particle Green's functions G°) identical to the equation 



B. The selfenergy: 

The vertex function can now be used to build up a full 
Green's function and different degrees of selfconsistency 
are possible. In this work we only want to focus on two 
ways to close the equations. 



1. Non self consistent calculation 

The first one is the non-selfconsistent calculation (see 
fig- §: 




FIG. 1. Diagram for the single-particle Green's function 
(solid line) in the non self-consistent approximation. The thin 
solid lines represent the non-interacting Green's function, the 
thick lines fully interacting Green's functions and the wavy 
lines represent the interaction U. 



G{k,iuj n ) = (G^k^)- 1 -E°(k,iw n ))" 



(5) 



were the selfenergy £°(k, iu> n ) is build up from the vertex 
function: 

£°(k,iw„) = — T° (k + q, iu m + ioj n )G° (q, iu m ) 



m,q 



(6) 



The system now consists of pairs of electrons which can 
only thermally be excited and can at finite temperatures 
be described as a mixture of Bosons (pairs) and Fermions. 
At zero temperature Bose condensation of the pairs into 
the two-particle bound state always occurs. This has 
been shown by Schmitt-Rinck et al. B were a version 
of the T-matrix theory was applied which was not even 
conserving on a one-particle level. 
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2. selfconsistent calculation 



The second possibility we want to discuss here is the 
fully selfconsistent version of the T-matrix. Other de- 
grees of selfconsistency have been discussed e.g. by Janko 
et al. J5J. Why we concentrate our discussion on the fully 
selfconsistent T-matrix calculation will be investigated 
elsewhere S. 




FIG. 2. Diagram for the single-particle Green's function in 
the self-consistent, conserving approximation, with the same 
notation as fig. 1. 

The full Green's function is now given by (see fig. ||): 

G(k,tw„) = (G (k,to n )- 1 -S(k,iw n )) _1 (7) 

were the selfenergy E(k, icu n ) is build up from the vertex 
function which is formed from a susceptibility being a 
product of full one-particle Green's functions as well: 



S(k, iu n ) 



—y 

N/3 ^ 



r(k + q,iu; m + iuj n )G(q,iu m ) (8) 



In the zero density limit (n/U = 0) eqs. (||) and (0) are 
identical. However at finite densities a different physical 
scenario is been described by the selfconsistent calcula- 
tion. This results from the fact that pair-pair interac- 
tions are now included in eq. (0). In this way the pairs 
can already at zero temperature have finite lifetime, they 
therefore deviate from Bose statistics and therefore the 
Bose condensation is strongly hindered (see fig. 0). 
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FIG. 3. The chemical potential is plotted as a function 
of temperature. The calculation was done for the interac- 
tion U = — St on a 12 by 12 lattice using 150 frequency 
points along the real axis. The dashed lines are results from 
non-selfconsistent calculations. At low temperatures they al- 
ways show Bose-condensation of pairs into the two-particle 
bound state below the line given by the non-selfconsistent 
Thouless criterion. In the selfconsistent calculation we al- 
ways regain a Fermi-like surface, but never reach a Thouless 
instability. For this plot the Hartree term was left out. 



C. Calculation procedure 

In order to solve the system of selfconsistency equa- 
tions we use a spectral representation for all correlation 
functions. The vertex function is turned into an analytic 
function T(K, iCl n ) by subtracting the Hubbard energy U 
which is identical to subtracting the Hartree term from 
the self energy. 



r(K,in n ) = f(K,in n ) - u 



u 2 x(K,in n ) 
i - u X (K, m n 



(9) 



We approximate the functions G(k, MJ n ), x(K, iQ n ), 
T(K,ifl n ) and S(k, iu> n ) by a scries of delta functions 
along the real axis e.g.: 



G(k,iu> n ) = 2J 



(10) 



to obtain a selfconsistent solution we have to determine 
the coefficients in an iterative way until selfconsistency 
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is achieved. The details of the method, especially how to 
choose the frequency points bj are explained in appendix 
A. 



IV. RESULTS 

We have solved the equations on a 2D system, a finite 
lattice with periodic boundary conditions. The results 
we present in the following were obtained on a 12x12 lat- 
tice. The attractive interaction was chosen to be equal 
to the bandwidth U = —8t. At finite densities we find 
completely different physics if we compare the selfconsis- 
tent and the non-selfconsistent calculations. In fig. ||| we 
have plotted the chemical potential as a function of Tem- 
perature. Note that the Hartree term is still ignored. In 
the non-selfconsistent calculation the chemical potential 
always ends up in the two-particle bound state, roughly 
-\y/U 2 + A 2 (with A the bandwidth) below the middle 
of the unperturbed band. If we solve the T-matrix self- 
consistently we get the chemical potential back into the 
one-particle continuum, we regain a Fermi-like surface 
but we never reach a superconducting instability. Note 
that the Thouless criterion in fig. || refers only to the su- 
perconducting instability of a non selfconsistent calcula- 
tion. Wc further get quasiparticle peaks at kp in the one- 
particle density of states whereas the non-selfconsistent 
calculation shows (at low enough temperatures) a gap at 
the chemical potential which is often interpreted as the 
origin of the pseudo-gap ffl]. For a large temperature 
range we further find a linear dependence of the quasi- 
particle scattering rate with temperature j9) which may 
explain the linear resistivity. 

In the current work we want to focus on two-particle 
properties. In the non-selfconsistent version of the T- 
matrix a two-particle bound state is present. At K = 
(tt, 7r) it evolves into which is called the ?7-mode p0| . This 
is shown in fig. [| (a). 









?7-mode: U • 


- 














1 „ " V 





FIG. 4. Schematic of the dispersion of the two-particle 
bound state below the non-interacting continuum of the sus- 
ceptibility x(K, Q.) (shaded region) as it arises in a non self- 
consistent formulation (a). At low enough temperatures the 
chemical potential ends always at the K = two-particle 
bound state. Below (b) we have plotted the same picture as 
it arises from a selfconsistent calculation done for the density 
n=0.4, temperature ksT=0.16[t] and interaction U = — St on 
a 12x12 lattice. The two particle bound state is disappeared 
and only a strong weakly dispersive two-particle resonance 
directly above the chemical potential (dashed line) is present. 

In order to answer the question how remnants of the 
two-particle bound state show up in the selfconsistent 
calculation we have plotted fig. [| (b). In this figure 
we show that the only remnant of the bound state is a 
weakly dispersive, strongly lifetime broadened resonance 
in the two-particle correlation function directly above 
the chemical potential. It is remarkable that also when 
Boson-Fermion models are discussed the Bosonic energy 
level is mostly put directly above the chemical poten- 
tial Also the 77-mode only partially survives the 
selfconsistent calculation. Even at K = (tt, tt) the two- 
particle resonance is strongly lifetime broadened and lies 
now less than U below the middle of the unperturbed 
band (here zero) . At zero density however we regain the 
ry-mode and the two solutions for the selfconsistent and 
non-selfconsistent calculation fall together. 
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FIG. 5. Lifetime of the pair for the selfconsistent calcula- 
tion above with the same parameters as in figs. 3 and 4. We 
have compared two different densities. With decreasing den- 
sity the lifetime of the pairs strongly increases until it reaches 
at n=0 infinity were also the 77-mode will be recovered. 

To quantify how strongly the pairs are lifetime- 
broadened we have plotted in fig. ^| the dispersion of 
the pair lifetime which is at a first approximation given 
by: 
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U 2 x"(O ) 



(11) 



were fio is the frequency at which the two-particle reso- 
nance occurs. The positions of the poles and the imag- 
inary part of the susceptibility at that points were ex- 
tracted from our numerical results. At K = (0, 0) the 
pairs live longest. With decreasing particle number the 
lifetime of the pairs increases until it reaches infinity at 
n=0 were the selfconsistent and the non-selfconsistent 
solutions fall together. 



V. CONCLUSION 

We have shown in this work that a selfconsistent in- 
clusion of pair-pair scattering into the T-matrix calcu- 
lation strongly alters the physics which is described by 
this theory. Due to selfconsistency the pairs can interact 
with each other which causes them to have finite lifetime. 
These finite lifetime pairs can no longer Bose condensate 
and therefore Fermi-liquid-like properties are regained. 
Although similar tendencies have been found earlier by 
other authors fl^|l4f| we are able to work out many de- 
tails and to reach very small temperatures. This is due 



to the fact that we use a semi-analytical method were all 
frequency integrations are done analytically and only the 
k-space summation has to be done numerically. Since we 
find the same results for finite lattices and for the con- 
tinuum in infinite dimensions Jl5| we do not believe that 
our results in general are sensitive to finite size effects. 
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APPENDIX A: SELFCONSISTENT PROCEDURE 

In order to solve the selfconsistent equations we apply 
a method which enables us to work entirely along the 
real frequency axis. We do an Ansatz for the correlation 
functions (here e.g. G(k, iuj n )) of the following form. 



G(k, iuj 71 



N 

E 



(12) 



This means the spectral representation of G(k, iu> n ) is 
approximated by a series of 5 functions. In this work 
the frequencies bj were kept fixed throughout the whole 
calculation (opposite to our previous work p"5| ) . In order 
to access low temperatures and work with a finite number 
of frequencies we sampled the frequency points with a 
tanh function. In this way we make sure that at the 
chemical potential two frequency points are always closer 
to each other then k^T. 

On the example of the susceptibility we want to show how 
we can calculate the spectral representation of a product 
function from the one-particle Green's functions: 



X (K, tfi n ) = _J_ £ G(K - k, iQ n 

m,k 



iu) m )G{k, iu) m ) 



(13) 

When inserting the spectral representation for G(k, iu) n ) 
we get: 

N 1 a K ~ k a k 

X<K E E j3 ^ iuj m - bj iQ, n - iuj m - bi 

k j,l m J 



N a*- k a k 



k j,l 



1 



1 



fin -bj-bi V 1 + e~f 3b i 1 + e~P b i 



9 X/X/ i 



N af - k a? 



2 ^ «fi n - bj - bi 



(14) 
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We now have determined x(K, iil n ) on N(N + l)/2 fre- 
quency points. These have to be sorted numerically and 
distributed onto the nearest frequency points bj . In order 
to calculate r(K,ir2„) from x(K, zO„) we have to apply 
a numerical broadening which must again be smallest at 
the chemical potential. In this way we get a complex 
function x(K,ii7„) from which we calculate the com- 
plex function r(K, iVL n ). From the imaginary part of 
r(K,ifi„) we calculate the amplitudes for the spectral 
representation for r(K,if2„). Finally we get: 

N K 

3 J 

A similar set of equations has to be solved for S(k, icu n ) 
and G(k, iu n ). In this way the selfconsistcncy is closed. 
We defined our solution to be selfconsistent if the mean- 
square deviation of all the amplitudes between two 
consecutive selfconsistency steps was below a certain 
threshold. 
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